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The scattering cross sections of charm mesons with hadrons such as the pion, 
rho meson, and nucleon are studied in an effective Lagrangian. In heavy 
O '. ion colhsions, rescattering of produced charm mesons by hadrons affects the 

', invariant mass spectra of both charm meson pairs and dileptons resulting 

Cf~) ' from their decays. These effects are estimated for heavy ion colhsions at SPS 

energies and are found to be significant. 
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> 
m 
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O ■ I. INTRODUCTION 

o . 

0^ ■ Recently, experiments on heavy ion collisions at CERN SPS by the HELIOS-3 and 

^ ! NA50 Collaborations have shown an enhanced production of dileptons of intermediate 
masses (1.5 < M < 2.5 GeV). In one explanation, this enhancement is attributed to dilepton 
production from secondary meson- meson interactions 0, while in another it was proposed 
that dileptons from charm meson decays could also contribute appreciably in this invariant 
Ch ■ mass region [Q. In the latter case, one of the present authors has shown, based on a 
schematic model, that if one assumes that the transverse mass spectra of charm mesons 
become hardened as a result of final state rescattering with hadrons, the invariant mass 



' spectrum of dimuons from decays of charm mesons would also become hardened, and more 
^ dimuons would then have an invariant mass between 1.5 and 2.5 GeV. Although charm 

quark production from hadronic interactions has been extensively studied using perturbative 
QCD [0,0, not much has been done in studying charm meson interactions with hadrons. 
Knowledge on charm meson interactions with hadrons is important as whether charm mesons 
develop a transverse flow depends on how strongly they interact with other hadrons as they 
propagate through the matter. 

In this paper, we shall first introduce in Sec. II an effective Lagrangian to describe the 
interactions of charm mesons with pion, rho, and nucleon. Using the coupling constants and 
cutoff parameters at the vertices determined either empirically or from symmetry arguments, 
we evaluate the scattering cross section of charm mesons with hadrons. Effects of hadronic 
scattering on the charm meson transverse momentum spectrum and the dimuon invariant 
mass spectrum from charm meson decays are then estimated in Sec. Ill based on a schematic 
model for the time evolution of heavy ion collision dynamics. In Sec. IV, we summarize our 
results and discuss the uncertainties involved in the studies. 
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II. CHARM MESON INTERACTIONS WITH HADRONS 



A. Lagrangian 



We consider the scattering of charm mesons {D~^, , D^, D^, D*^, D*~~ , and 
D*^) with pion, rho, and nucleon. If SU(4) symmetry were exact, interactions between 
pseudoscalar and vector mesons could be described by the Lagrangian 

Cppv = igTr [P^V^d^P) + H.c. (1) 

where P and V represent, respectively, the 4x4 pseudoscalar and vector meson matrices 
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The above interaction Lagrangian may be considered as being motivated by the hidden 
gauge theory, in which there are no four-point vertices that involve two pseudoscalar mesons 
and two vector mesons. This is in contrast to the approach of using the minimal substitution 
to introduce vector mesons as gauge particles, where such four-point vertices do appear. It 
is, however, known that the two methods are consistent if one also includes in the latter 
approach the axial vector mesons, which are unfortunately not known for charm hadrons. 
Furthermore, gauge invariance in the latter approach cannot be consistently maintained if 
one uses the experimental vector meson masses, empirical meson coupling constants and 
form factors at interacting vertices. Expanding the Lagrangian in Eq. (0) in terms of the 
meson fields explicitly, we obtain the following Lagrangians for meson-meson interactions: 



^ttDD 
^pDD 

C 



'- WttDD 

WpDD 
7^ 
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H.c. 



Drid^D) - {d^D)rD 



(2) 



where the coupling constants g-jrDD*, dpDD, and (/p^^r are related to the coupling constant g 
via the SU(4) symmetry as shown below in Eq. (^. 

In SU(3), the Lagrangian for meson-baryon interactions can be similarly written using 
the meson and baryon matrices. The formulation becomes, however, more complicated in 
SU(4) where a more general tensor method is required 0. The interaction Lagrangians 
needed for our study then include the following: 

C^NN = -ig^NNN-i^rN ■ 7? , 

C-DNK = Wdnk, {n^^K^D + Ac75A^^) • 
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In the above we have used the following conventions: 

TT^ = ^^1^, N = {p, n), and D = (D^, 

Again, SU(4) symmetry would relate the above coupling constants to each other with the 
introduction of one more parameter, as shown below in Eq. (^, because there are two SU(4)- 
invariant Lagrangians for pseudoscalar meson and baryon interactions. We also need the 
following phenomenological Lagrangian: 

'C.pNN = QpNnN (^7^f ■ + ^^o-'^^r ■ d^p^ N , 

where values of the coupling constants QpNN and Hip are well-known as discussed below. 



B. Cross sections 



In Fig. |I|, Feynman diagrams are shown for charm meson interactions with the pion 
(diagrams 1 to 8), the rho meson (diagrams 9 and 10), and the nucleon (diagrams 11 to 13). 
Explicit isospin states are not indicated. The spin and isospin-averaged differential cross 
sections for the t channel and u channel processes can be straightforwardly evaluated, and 
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where Pi denotes the initial momentum of the two scattering particles in their center-of-mass 
frame. 

For s channel processes through charm meson resonances, shown by diagrams 4 to 8, the 
cross section is taken to have a Breit-Wigner form 



(2 J + 1) AtT Ft^,, Angput 

(2.i + l)(2s2 + l) pUs-Ml^y/s + Tl 



3 



where Ftot is the total width of the resonance, -Bin and -Bout are their decay branching ratios 
to the initial and final states, respectively. We note that diagrams 4 to 7 correspond to 
processes through the -Dg and Di resonances Dir -Dtt, -Dtt — -D*7r, D*tt Dtt and 
D*7r D*7r, respectively, while diagram 8 represents the process Dtt — ^ Dtt through the 
D* resonance. Total widths for -Dg and Di resonances are known, and they are r£)*o = 23 
MeV, r^.+ = 25 MeV, T^o = 18.9 MeV, and = 28 MeV g. For the width of D*, only 
an upper limit is known, i.e., F^i.o < 2.1 MeV and T£)*+ < 0.131 MeV. Studies based on the 
relativistic potential model suggest that T£)*o ~ 42 KeV and T£i*+ ~ 46 KeV, and we 
use these values in this paper. The branching ratios (BR) are known for D*^ and D*^ [Q, 
but not for -Dg and Di. Experimental data show that for both -Dg^ and D'^'^ decays one has 
T{DTTch) /^{D*7Tch) ~ 2. Since Di decays to D*tt instead to Dtt due to parity conservation, 
we assume BR{Di D*tt) = 1, BR{D* D*tt) = 1/3, and BR{D* Dtt) = 2/3, 
neglecting possible decays of Di and -Dg to Dp and D*p, respectively [p!0|j . 



C. Coupling constants 

For coupling constants, we use the empirical values gp-^^ = 6.1 [|n[], g-^DD* = 4.4, gpon 
2.8 |T^, g-j^NN = 13.5 |1^, gpNN = 3.25, and Kp = 6.1 H^. From SU(4) symmetry, as 



assumed in the Lagrangian in Eq. (P, one would expect the following relations among these 
couplings constants: 

9.KK*i^-^) = ^7.DD*(4.4) = gpKKim = 9pDd{2.8) = ^(3.0) . (4) 

One sees that the empirical values given in parentheses agree reasonably well with the 
prediction from SU(4) symmetry. Signs of the coupling constants are not specified as the 
possible interferences among diagrams 3, 4, and 8 are not included. We note that the 
coupling constant g^oD* is consistent with that determined from the D* width 

■p _ dloD'-P} 

J- D*-*TrD — — 2 ' 

2TTmf), 

where pj is the momentum of final particles in the D* rest frame. 



D. Form factors 



To take into account the structure of hadrons, we introduce form factors at the vertices. 
For t channel vertices, monopole form factors are used, i.e.. 



m 



A^-t ' 

where A is a cutoff parameter, and rua is the mass of the exchanged meson 
parameters, we use the empirical values Ap^rTr 
1.4 GeV 



1.6 GeV Q, A 



tNN 



For cutoff 
1.3 GeV, and ApTVA^ = 



T3|. However, there are no experimental information on K^^dd* and Ap^ij, and their 
values are assumed to be similar to those determined empirically for strange mesons, i.e., 
AttDD* = At,kk* = 1-8 GeV, ApCD = ApXK = 1-8 GeV [0. For s channel processes, 
shown in diagrams 4 to 8, that are described by Breit-Wigner formula, no form factors are 
introduced. 
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E. On-shell divergence 



The cross sections in Eq. @ for diagrams 2 and 9 {Dp ^ D*7i) are singular because the 
exchanged mesons can be on-shell. Since the on-shell process describes a two-step process, 
their contribution needs to be subtracted from the cross section. This can be achieved by 
taking into account the medium effects which add an imaginary self-energy to the mass of 
the exchanged pion as in Ref. [15]. We take the imaginary pion self-energy to be 50 MeV 
and have checked that the calculated thermal average of the cross sections do not change 
much with values between 5 and 500 MeV. We note that there are other ways to regulate 
this singularity ||16 . 



F. Thermal average 

We are interested in the thermal averaged cross sections for the processes shown in Fig. |l|. 
For a process 1 + 2 — >■ 3 + 4, where the initial-state particles 1 and 2 are both described by 
thermal distributions at temperature T, the thermal averaged cross section is given by 

. . j:;^dz[z'-{ai + a2f][z'-{ai-a2f]K,{z)a{s = z'T') 

(crv) = — . 

4(1 + 5i2)aii^2(«i)a2^2(a2) 

In the above, = rrii/T (z = 1 to 4), zq = max(ai + ^2, Q^s + "4), ^12 is 1 for identical 
initial-state particles and otherwise, and v is their relative velocity in the collinear frame, 
i.e., 

V (ki ■ — mfml 

V = . 

E1E2 

In Fig. H(a), we show the results for the thermal averaged cross sections as functions 
of temperature. It is seen that dominant contributions are from D and D* scatterings by 
nucleon, D scattering by pion via rho exchange, and D scattering by rho meson via pion 
exchange. In obtaining these result, the rho meson mass is taken at its peak value of 770 
MeV. 



III. ESTIMATES OF RESCATTERING EFFECTS 

As shown in the schematic model of Ref. [Q, if one assumes that charm mesons interact 
strongly in the final state hadronic matter, then their transverse mass (m^) spectra and 
pair invariant mass spectra would become harder than the initial ones as a result of the 
appreciable transverse flow of the hadronic matter. Dilepton decays of charm mesons would 
then lead to an enhanced yield of intermediate-mass dileptons in heavy ion collisions. In 
this section, we estimate the effects of hadronic rescattering on charm meson spectra 
and the invariant-mass distribution of dileptons from their decays in heavy ion collisions at 
SPS energies. 

To characterize the scattering effects on charm mesons, we first determine the squared 
momentum transfer to a charm meson when it undergoes a scattering process DiXi D2X2. 
In the rest frame of Di, the squared momentum of the final charm meson D2 is given by 
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for t channel processes with four momentum transfer t. For u channel processes, one replaces 
u for t in the above expression. 

We determine the total number of collisions suffered by a charm meson from its scattering 
cross sections and the time evolution of the hadron densities. In the charm meson local frame, 
we assume that the density evolution of hadrons is inversely proportional to the proper time, 
i.e., 

p(r) oc - . (5) 

r 

Neglecting the effect of transverse expansion on the density evolution, the total number of 
scatterings for a charm meson is then 

f^F ( Tp 

N = avpdr = avp^T^ In — 
Jto Vto 

= avpoTo In — ~ avpoTo In , 

\TocoshyJ \ Top± J 

which leads to the following thermal average of the squared total momentum transfer due 
to scatterings. 



{pD = {NpD 



{(^'"pDipio 

=n,p,N--- 



r„ln(^). (6) 



In obtaining the above result, we have assumed the same initial and final proper times for the 
time evolution of different particle species that are involved in the scattering. Equation 
shows that the relevant quantity is (crvpl) instead of the usual (erf). We show in Fig. ^(b) 
this thermal average for all scattering channels considered in the present study. It is seen that 
the dominant contributions to (av) remain important for (crvpl) and the process involving 
D* scattering by rho meson via pion exchange also becomes significant. 

Summing up contributions from the scattering channels in Fig. [^(a), |l](b), and |l|(c) 
separately, and simply dividing by 2 to account for the average over D and D*, we obtain, 
at T = 150 MeV, 

(avpl) - 1.1, 1.5, and 2.7 mb GeV^ 

for vr, p, and scatterings with charm mesons, respectively. 

For central Pb+Pb collisions at SPS energies, the initial particle numbers can be obtained 



from Ref. 0, i.e., there are 500 vr, 220 p, 100 80 r], 180 N, 60 A, and 130 higher 
baryon resonances. The initial densities at central rapidity can then be estimated using 
PqTq ~ {dN/dy)/{7iR\) ~ N/{4:nR\). For a conservative estimate on the scattering effect, 
we only include tt, p, and nucleon. The initial densities for pion, rho meson, and nucleon 
are thus 

PqTq ~ 0.79, 0.35, and 0.28 fm'^ , 
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respectively. Equation (D then gives 



~ (1.1 X 0.79 + 1.5 X 0.35 + 2.7 x 0.28)/10 x In 16.7 - 0.61 (GeV^) . (7) 
In the above, we have taken tq = 1 fm and R± ~ Ra — 1.2A^^^ fm. We have also used 



the relations {p±) ~ y {p\) ^ v2m7^ and (m^) ~ m + T^e as given by Eq. (^) in 
Appendix A. Since the charm meson Teg increases as a result of the rescatterings, (pg) needs 
to be determined self-consistently. However, because of the logarithmic dependence shown 
in Eq. (^), (pg) is not very sensitive to the value of Tefj, and we have taken Teg = 200 MeV 
in obtaining the above numerical results. We note that even though pions appear to be less 
important in Fig. ^ their contribution to the rescattering effect is important due to their 
high densities, as evident from the numerical values shown in Eq. (0). 

The total squared momentum transfer from hadronic scatterings as given by Eq. can 
be characterized by a temperature parameter Tg, defined by Eq. ( |A^ ) in Appendix A. Using 
the values given in Eq. (|^), we obtain Tg ^ 96 MeV from Eq. ([A5|) of Appendix A. From 
Fig. 1^, which relates Tg to Teg and to the enhancement factor R for dimuons from charm 
meson decays into the NA50 acceptance, this gives an effective inverse slope parameter of 
Tgff = 235 MeV for the final charm meson m±_ spectrum if the initial one is taken to be 160 
MeV, and a dimuon enhancement factor of about 2.1 is obtained. 



IV. SUMMARY AND DISCUSSIONS 

In summary, we have calculated the cross sections for scatterings between charm mesons 
and hadrons such as pion, rho meson, and nucleon. Hadronic scatterings of charm mesons 
in heavy ion collisions can significantly affect the charm meson spectra and the dilepton 
spectra from charm meson decays. An estimate of this effect in heavy ion collisions at SPS 
energies is given, and we find that it leads to a hardened charm meson spectra and an 
enhanced intermediate-mass dileptons from charm meson decays. These results thus give a 
more quantitative justification of the arguments proposed in Ref. [Q. 

However, the results obtained in the present study are still incomplete as we have not 
included diagrams involving the exchange of heavier particles such as charm hadrons. The 
scattering cross sections between charm mesons and hadrons such as kaon, u;, r/. A, and 
higher baryon resonances are not calculated either. Furthermore, we have not calculated the 
contribution due to diagrams shown in Fig. ^, where charm mesons scatter with nucleons via 
a Ac exchange. There is a large uncertainty in their contributions as no empirical information 
on the coupling constant Qdnk^ is available. The SU(4) symmetry gives 

3 — 2a/) , . 

QDNAc = 9KNk = QirNN — dnNN = 13.5 , (8) 

where an = D/{D + F) with D and F being the coefficients for the D-type and T-type 
coupling, and ~ 0.64 ||18|. On the other hand, QCD sum rule studies suggest a smaller 
value QdnAc — 6.7 ± 2.1 W^. Because of this uncertainty in these two processes involving 
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the DNAc coupling, we choose to leave them out in our study. In the future, when empirical 
information on this parameter is known from DN scattering, then processes involving the 
DNAc coupling can be addressed. 

We note that the estimates given above are based on a simple assumption on the time 
evolution of the dense hadronic system, which enables us to make an analytical estimate 
of the rescattering effects. As a result, we have neglected the transverse expansion of the 
hadronic system which would lead to a faster decrease of hadron densities than the linear 
dropping assumption in Eq. (|]). We have also neglected the chemical equilibration processes 
which, e.g., may decrease the total number of rho mesons and increase that of pions as a 
function of time |T7 |. 

Moreover, we have used only the isospin averaged cross sections and also averaged the 
rescattering effects on D and D* mesons. Without a full cascade calculation and fully 
treating the isospin, we do not know the final composition of charm mesons, e.g, the ratios 
D*/D and D^/D~^. A naive expectation gives D*/D = 3, and consequently D^/D^ ~ 3 
20| . However, even for pp collisions the relative weights of produced charm mesons are not 



well measured experimentally. We emphasize that the charm meson composition could have 
a sizable effect on the lepton and dilepton yields from charm decays, because and 
have very different branching ratios for semileptonic decays (17.2% from and 6.7% from 
DO). 

In a hadronic cascade model, the time evolution and the chemical equilibration of the 
hadronic system can be simulated much better. Using cross sections with the full isospin 
information, and keeping track of the charm meson isospins during scatterings, the final 
charm meson composition can be determined. Therefore, further studies based a cascade 
code along these directions are much needed for a quantitative study of the rescattering 
effects on charm meson observables. 

For heavy ion collisions at RHIC energies, a dense partonic system is expected to be 
formed during the early stage of the collision. In addition to hadronic rescatterings of charm 
mesons, partonic rescattering effects on charm quarks also need to be included. Furthermore, 
radiative processes of charm quarks inside the QGP would further complicate the issue as 
they may cause energy loss and soften the charm meson m± and pair invariant mass 



spectra [22]. Therefore, more studies are needed before one can make predictions for RHIC. 
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APPENDIX 

In this appendix, we derive the relation between the total squared momentum trans- 
fer to charm mesons due to hadronic scatterings and the increase of the inverse slope of 
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charm meson m± spectra. Consider a charm meson at central rapidity with an initial trans- 
verse momentum p_i_i along the y axis. After a scattering which gives the charm meson a 
momentum ps in its rest frame, its final transverse momentum is given by 

PxF = PxS , PyF = lliPyS + A^s) , 

where 

a P±l 
A = ; , , 7l 



Assuming that ps is isotropic in the charm meson rest frame, the average of the squared 
final transverse momentum of the charm meson is then related to that of the squared initial 
transverse momentum by 



(Pi.) = (Pi.) + (| + ^) (pD , (Al) 



where (pg) is the average of the squared total momentum transfer to the charm meson as 
given by Eq. (^. For an isotropic ps distribution, Eq. ( |Al| ) is actually true for a charm 
meson at any rapidity. 

If we parametrize the spectrum of charm mesons as 



m^dm 

in terms of an inverse slope parameter Teg, then 

{pD = 2rJr f + 2 + -j.^) ■ (A2) 
\Tefr m/Tcff + iy 

As in the schematic study of Ref. we characterize the scattering strength (pg) by an 
equivalent temperature parameter Ts via 



2, Ip'e-^/^^d^p ^2 /^^^P (-V^' + {m/Tsf)dx 



(A3) 



3T| 



m 5 / 1 
— — — 

Ts 2 \m/Ts 



Both Ts and T^s are expected to be small compared with the charm meson mass (m ~ 
1.87 GeV). Keeping only the leading term in Eqs. ( [A2| ) and ( [A^ ) then gives = 2mTeff 
and (pg) = 3mTs in the nonrelativistic limit. Equation ( [Al] ) thus gives 

Tjff ^ Tifj + Ts . (A4) 

If we also keep the next-to- leading term in Eqs. ( |A2D and (|A3|), we then obtain 

{pD ^ 3mTs + ^T| (A5) 
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For central Pb+Pb collisions at SPS energies, the relations given by Eqs. ( |A4| ) and ( |A6| ) 
are shown in Fig. ^ together with the results obtained from the Monte Carlo simulations in 
the schematic study of Ref . @] . These relations agree qualitatively with that from the sim- 
ulations. However, they differ quantitatively, because rapidity changes due to rescatterings 
are not taken into account in the present analytical estimates. 



10 



REFERENCES 



M. Masera for the HELIOS-3 Collaboration, Nucl. Phys. A590, 93c (1995). 
E. Scomparin for the NA50 Collaboration, Nucl. Phys. A610, 331c (1996). 
G. Q. Li and C. Gale, Nucl. Phys. A638, 491c (1998). 
Z. Lin and X.-N. Wang, Phys. Lett. B 444, 245 (1998). 

B. L. Combridge, Nucl. Phys. B151, 429 (1979). 

P. Nason, S. Dawson, and R. K. Ellis, Nucl. Phys. B303, 607 (1988); B327, 49 (1989); 
W. Beenakker, H. Kuijf, W. L. van Neerven, and J. Smith, Phys. Rev. D 40, 54 (1989). 
S.K. Bose and W.D. McGlinn, Phys. Rev. D 14, 3167 (1976). 
Review of Particle Physics, Particle Data Group, Euro. Phys. J. C 3, 1 (1998). 
P. Colangelo, F. De Fazio, and G. Nardulli, Phys. Lett. B 334, 175 (1994). 

E. J. Eichten, C. T. Hill, and C. Quigg, Phys. Rev. Lett. 71, 4116 (1993); Y.-B. Dai 
and H.-Y. Jin, Phys. Rev. D 52, 236 (1995). 

C. M. Ko and D. Seibert, Phys. Rev. C 49, 2198 (1994); G. E. Brown, C. M. Ko, Z. G. 
Wu and L. H. Xia, Phys. Rev. C 43, 1881 (1991). 

S.G. Matinyan and B. Miiller, Phys. Rev. C 58, 2994 (1998). Note that our coupling 
constant is a factor of 2 smaller due to the difference of the definitions. 
B. Holzenkamp, K. Holinde, and J. Speth, Nucl. Phys. A500, 485 (1989); G. Janssen, 
J. W. Durso, K. Holinde, B. C. Pearce, and J. Speth, Phys. Rev. Lett. 71, 1978 (1993). 
G. Janssen, K. Holinde, and J. Speth, Phys. Rev. C 54, 2218 (1996). 
G. Q. Li and C. M. Ko, Nucl. Phys. A594, 439 (1995); W. Peters, U. Mosel, and A. 
Engel, Z. Phys. A 353, 333 (1996). 

R. Baier, M. Dirks, and K. Redhch, Phys. Rev. D 55, 4344 (1997). 
G. Q. Li, C. M. Ko, G. E. Brown, and H. Sorge, Nucl. Phys. A611, 539 (1996). 
R. A. Adelseck and B. Saghai, Phys. Rev. C 42, 108 (1990). 

F. S. Navarra and M. Nielsen, Phys. Lett. B 443, 285 (1998). 
S. Frixione, M. L. Mangano, P. Nason, and G. Ridolfi, |hep-ph/9702^ . 
R. Baier, Yu. L. Dokshitser, A. H. Mueller, S. Peigne, and D. Schiff, Nucl. Phys. B478, 
577 (1996); B483, 291 (1997); R. Baier, Yu. L. Dokshitzer, A. H. Mueller, and D. 
Schiff, Phys. Rev. C 58, 1706 (1998). 

[22] Z. Lin, R. Vogt, and X.-N. Wang, Phys. Rev. C 57, 899 (1998). 



11 




12 




4^ 




100 120 140 160 180 



120 140 160 180 200 



T (MeV) 



FIG. 2. Thermal average (a) {av), and (b) {(tvpq) of charm meson scattering cross sections 
as functions of temperature. Numbers labeling the curves correspond to the diagram numbers in 
Fig-S 
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FIG. 3. The dimuon enhancement factor R within the simulated NA50 acceptance and the 
equivalent temperature parameter To, due to scatterings as functions of the final inverse slope of 
charm mesons Tefj Q. The curve with open circles is from Eq. (X4), while the curve with filled 
circles is from Eq. (|A6|). 
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